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A STUDY OF HIGH ORDER VERTICAL APERTURE CORRECTION 
J.O. Drewery, M.A., Ph.D. 



1. Introduction 

The loss in vertical resolution of television cameras and 
teleclne machines caused by finite spot size is commonlv 
compensated for by using a vertical aperture corrector. 
This operates by modifying the signal associated with any 
particular line using information derived from signals on 
adjacent lines and, until now, the correction information 
has been derived from the signals on only the two nearest 
neighbouring lines. The amount of correcting signal can be 
varied to suit the particular source under consideration but 
if too much correction is applied, horizontal edges in the 
displayed picture are accompanied by overshoots and, in 
extreme cases, by 'ringing' extending for many lines above 
and below the edges. 

The object of the work presented in this report was 
to investigate the effect of applying more complex correc- 
tion using information derived from many adjacent lines in 
the hope that by so doing an improved degree of subjective 
sharpness could be achieved without the impairment 
caused by over^oot. 

The penalty of increased complexity, apart from cost, 
is the increase in noise introduced into the corrected signal 
as more and more components contribute to the final 
result However, this is less important than it would seem 
becajse the additional noise components are weighted by 
coefficients wh\c\n fall off very rapidly as the complexity 
increases. This will be shown below. 

The application of more complex correction requires 
a means of simultaneous access to the signals on several 
adjacent lines. This may be achieved by using an arrange- 
ment of cascaded line delays. Analogue delays of this 
length always produce small but measurable signal impair- 
ments and such impairments have, so far, limited the com- 
plexity of correction. Digital delays, however, provide a 
means of delaying signals with in principle no degradation 
and thus of realising correction of any degree of com- 
plexity. 

2 3 

The prototype digital line-store standards converter ' 
offered a means of experimenting with complex aperture 
correction, as very little modification to the equipment was 
required to make it operate as an aperture corrector. The 
modification consisted simply of taking the output directly 
from the interpolator and adding digital signals correspond- 
ing to the synchronising and blanking signals before feeding 
the signal to the d.a.c. The interpolator coefficient 
generator had also to be modified slightly so as to give a 
constant set of values obtaining throughout the whole 
field. 

Such an arrangement gives the experimenter many 
degrees of freedom and, if undue experimental time is not 
to be wasted, a basis of theory is required which can guide 
investigations in the dir«;tion of the parameter values most 



likely to achieve optimum results. Accordingly space is 
devoted here to a theoretical treatment of aperture correc- 
tion before practical experiments are described. 



2. Tlieoretical considerations 

2.1. The theory of aperture degradation 

A television system scans an object field with a finite 
aperture and in so doing converts a two-dimensional spatial 
function into a one-dimensional temporal function or signal. 
This signal is used at the display to modulate the light 
passing through a second aperture scanning in sympathy, so 
recreating approximately the original scene. The theory of 
scanning involving the use of finite apertures was first 
deduced by Mertz and Gray"* and is so central to the 
particular problem of aperture correction that it will be 
reviewed here. 

The object field can be considered as the superposition 
of two-dimensional waves having harmonically related 
spatial frequencies, much in the same way as a one- 
dimensional function can be broken down into its Fourier 
components. A general physical two-dimensional spatial 
frequency has, say, m cycles per picture width and n cycles 
per picture height and can be expressed as the sum of two 
complex spatial frequencies having indices m, n and -m, 
-n. There is however a further physical two-dimensional 
spatial frequency having m cycles per picture width and n 
cycles per picture height which has a slope of opposite sign. 
This can be expressed as the sum of two complex fre- 
quencies of indices m, -n and -m, n. A general complex 
frequency can be denoted simply by the pair (m, n) having 
an amplitude ^^„. The scanning action converts each 
complex spatial frequency (m, n) into a conventional tem- 
poral frequency, having an amplitude proportional to^^„ 
and a frequency /given by: 



/= ot/h + n/v 



(1) 



where /^ and /^ are the horizontal and vertical scanning 
frequencies respectively. If the picture is horizontally 
scanned sequentially in iV lines /^ =/l, the line frequency, 



f\i ^ fp' ^'^^ picture frequency, and f^_ 
Equation (1) becomes 

/=4(m + n/7V) 



NL. Then 



(2) 



If N is very large. Equation (2) shows that the signal 
spectral energy appears in characteristic bunches centred on 
multiples of the line frequency as shown in Fig. 1. 
Equation (2) also shows that although any given spatial 
frequency (m, n) transforms into one temporal frequency 
the reverse is not true. For any pair of values (m , n ) 
which satisfy the relationship 



m + n'lN =m + n/N 



(PH-118) 



0,0 



-1,0 



2,0 



3,0 



spatial frequency c/pw,c/ph 



LJ 



ZfL 



II Illlll 



3^1 



Ulu, 



temporal frequency, Hz 



Fig. 1 - Characteristic form of sequential television signal 
^ectrum 

also corresponds to the same temporal frequency. In 
other words, if the values of m and n are unrestricted the 
bunches of energy can overlap. This effect was termed 
'confusion' by Mertz and Gray. Nowadays this would be 
explained by saying that the scanning action spatially 
samples the scene in the vertical direction so that 'aliasing' 
can occur if the a;ene contains components having vertical 
frequencies greater than half the sampling frequency. Thus 
the aliasing or confusion can be avoided by restricting the 
vertical frequencies of the original scene to values of n 
below Nl2. This brings out an important distinction in the 
status of the horizontal and vertical directions, in that it is 
not necessary to restrict the values of m to avoid confusion. 
In other words there is no basic sampling in the horizontal 
direction. 

The effect of using a finite scanning aperture is to 
modify the amplitudes of the spatial frequency components 
in the scanned signal, related to A^^, by a factor Y(m,n) 
which is the two-dimensional Fourier transform of the 
reversed scanning aperture.* This is important because it 
means that it is possible to consider the object field to be 
pre-distorted by the aperture before scanning. Thus, for 
example, it is possible to avoid confusion in the signal by 
choosing a suitable scanning aperture which attenuates 
^atial frequencies having values of n greater than Nil. 
This is entirely analogous to the concept of appropriately 
restricting the baseband frequency range to avoid aliasing in 
a sampling process. 

If the original scene contains no components with 
values of n greater than NI2 then clearly the function of 
aperture correction is to restore the amplitudes of the com- 
ponents in the scanned signal to the values that would have 
been obtained had the scanning been carried out with a 
point aperture. If, however, the possibility of confusion 
exists it is a difficult matter to decide whether or not to 
restore the components in their entirety and, if not, whether 
to apply a sharp cut or gentle 'roll-off. 

In principle it is only necessary to pass the signal 
through a network whose response/frequency characteristic 
is the reciprocal of Y(m,n), n<Nl2, expressed in terms of 
temporal frequency. In practice this is done in two stages 
known as horizontal and vertical aperture correction. These 
two processes can, however, provide only an approximation 
to what is required in certain cases as will now be shown 



If the scanning aperture is T[x,y) the reversed aperture is T(—x,~y). 
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Fig. 2- Spectral transfer function of aperture loss 



The general effect of aperture degradation is to 
attenuate progressively the amplitudes of the higher fre- 
quency components. Thus a horizontal aperture corrector 
has a response which rises steadily throughout the video 
bandwidth and can be arranged to correct perfectly the 
components at multiples of line frequency i.e. the hori- 
zontal spatial frequencies (m,0). The form of vertical 
aperture corrector commonly used is the symmetrical trans- 
versal filter whose response/frequency characteristic is 
periodic at multiples of line frequency; thus if the corrector 
response is unity at zero frequency the spatial frequencies 
(m,0) are also uncorrected. Roughly speaking therefore 
the horizontal corrector affects the size or envelope of the 
spectral energy bunches and the vertical corrector the 
shape, as shown in Fig. 2. Now implicit in this approach is 
the assumption that the degradation of the shape of the 
bunches is uniform throughout the video bandwidth. This 
in turn implies that the aperture spectrum is variables- 
separable so that it can be expressed in the form:— 

Y(m,n)=YJm)Yy(n) 

where Y^(m) behaves as an envelope function and Yy{n) as 
a shape function independent of m. This assumption holds 
for the rectangular and elliptical Gaussian apertures, but 
not, for example, for the circular aperture. Thus, the con- 
ventional aperture correction process cannot correct per- 
fectly for a circular aperture or any aperture that cannot be 
expressed in variables-separable form. 

In order to correct for such apertures a generalised 
form of transversal filter would be required whose delay 
elements would consist of alternating groups of element 
delays and shortened line delays. Such an arrangement 
gives independent access to each picture point of a cluster. 

In the remainder of this discussion attention will be 
directed towards the vertical aperture correction process, 
assuming that it is justifiable to separate the horizontal and 
vertical effects. Indeed, it is only in thiscontext that the 
terms horizontal and vertical aperture have any real mean- 
ing. The process of vertical aperture distortion and correc- 
tion can then be treated as a one-dimensional spatial 
problem involving spatial sampling. 

2.2. The vertical aperture correction process 

The spatial field down a vertical line in the scene can 
be expressed as a superposition of complex sinusoids of 
amplitude A^ where the >l„'s are functions of the position 
of the vertical line and the amplitudes yl^„, that is 
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Fig. 3 - The correction process in trte spectral domain 

(a) original signal (6) source aperture (c) sampled signal 

id) correction aperture 

(e) convolution of (c) and id) in spatial domain 

if) corrected and resampled signal 



E{x^,y)= / .A„. exp]my/b 



where 



A 



„= /] A^„.eKp\7rnvcJa 



and the scene has dimensions 2a and 2b along the x andy 
directions respectively. After scanning with a finite 
variables-separable aperture the pre-distorted field down the 
same vertical line could be expressed as 






r (n) . exp]my/b 



where 



A^„ . Y^{m) . expjrrmXj/a 



Assuming that the horizontal degradation is subsequently 
corrected so that y4„' becomes ^„, the effect of vertical 
degradation is to multiply the vertical spectrum by the 
function Yy («), the spatial Fourier transform of the vertical 
aperture. Thus it is possible to frame the discussion in 
terms of the one-dimensional spectral functions A„ and 
Yy{n). 

In order to correct the degradation it is necessary to 
operate on the signal, which gives values of the degraded 
image at equispaced points down the vertical line, by 
restoring these values to those that would have obtained 
had the scanning been carried out with a point aperture. 
Remembering that sampling action results in a spectrum 
with sidebands centred on multiples of the sampling fre- 
quency, the sequence of operations involved in correction 
could be described in the spectral domain as in Fig. 3. The 
original signal whose vertical spectrum, ^„, is shown at (a) 
is scanned with an imperfect source aperture whose vertical 
spectrum, Y(n), is shown at (&) to produce a signal whose 
spectrum is shown at (c). Multiplication by the function 
shown at {d) produces the spectrum shown at (e) which can 
then be rescanned in the space domain to produce a signal 
approximating to perfect scanning whose spectrum is shown 
at if). 

It will be observed that the correction function shown 
at {d) departs from the reciprocal of the source aperture 
spectmm at frequencies above about half the sampling fre- 
quency. This must be so, for two reasons; first, as the 
source aperture spectrum falls to zero the reciprocal cannot 
tend to infinity and, secondly, the higher spectral orders 
must be eliminated if they are not to intrude into the zero- 
order spectrum after resampling (i.e. by aliasing). 

The correction process represented by the step from 
(c) to {/■) has a transfer function which is periodic in the 
spatial frequency domain and is, in fact, obtained by 
superposing versions of the function shown at (d) repeated 
at intervals of A^. Now the symmetrical transversal filter 
shown in Fig. 4(fl) can, with certain values of the coef- 
ficients, c^, be given the response shown in Fig. Mb). If 
the value of r is a line period the transformation of spatial 
to temporal frequency implied by Equation (2) means that 
there is clearly a relationship between the function shown 
at Fig. 3(d) and the coefficients c^, as will now be shown. 

The multiplication of the spectrum of Fig. 3(c) 
representing the degraded samples by the function shown at 
id) is equivalent to convolving the samples with a function 
which is the inverse Fourier transform of that shown at {d). 
That is to say, each sample is replaced by a scaled version of 
this inverse Fourier transform which, by analogy, could be 
called a correction aperture. The continuous function so 
produced is then resampled at the same sample points to 
yield modified samples. This process is shown in Fig. 5 
which is the spatial equivalent of Fig. 3. As, in general, the 
correction aperture will not have zeros at the sampling 
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Fig. 4 - The transversal filter 
ia) network (6) a possible transfer characteristic 



points, the magnitude of each sample will modify its 
neighbours. Now returning to the transversal filter of Fig. 
4{a), if the input signal is aich that the signal magnitude is 
zero at all the tapping points except one, as the signal 
advances steadily through the filter the output of the filter 
is that signal magnitude multiplied by successive coefficient 
values. But this is exactly what would have been obtained 
by convolving the isolated sample with the correction 
aperture and resampling the result at the sample points. 
Thus the coefficient values are the values of the correction 
aperture at the sample points. 

The corrector is adjusted by varying the values of the 
coefficients until the response over the zero-order spectral 
region approximates to [Yy («)] ~ ' . It will now be realised 
that in order to correct for an aperture perfectly up to the 
theoretical limit of displayed spatial frequency, that \sN/2 
cycles/picture height, the corrector response must have a 
cusp at this frequency because the response is periodic. 
This in turn, implies that there must be an infinite number 
of coefficients. Looked at another way, for perfect 
correction the function shown at Fig. 3{d) must fall dis- 
continuously to zero at the limit. Thus the correction 
aperture, being the inverse Fourier transform, will 'ring' 
indefinitely producing an infinite number of coefficients. 
Thus, not unexpectedly, a perfect aperture correction 
requires a corrector of infinite extent. 

Although the goal of perf«;t correction can never be 
achieved, it is nevertheless instructive to calculate the ideal 
coefficient values for an infinite corrector because they 
indicate the order of magnitude of the quantities involved. 



Fig. S - The correction process in the spatial domain 

(a) original signal ib) source aperture (c) sampled signal 

id) correction aperture (e) corrected signal 

(f) corrected and resampled signal 

2.3. Perfect smirce correction 

The preceding section indicates a method of deducing 
coefficient values. The steps are:— 

Find the spatial Fourier transform of the vertical 
aperture. 

Take the reciprocal and truncate it so that it extends 
over a region of width equal to the spatial Planning 
frequency. 

Take the inverse Fourier transform. 

Calculate values of the function at the sampling- 
(scanning) points. 

This method can be applied to any vertical aperture, 
provided that its spatial Fourier transform does not have a 
zero within the range ±N/2 c/p.h. Fig. 6 shows the values 
of the first few coefficients when the method is applied to 
rectangular and Gaussian apertures. The derivation of 
these curves is given in Appendix 1 where it will be seen 
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Fig. 6 ■ Theoretical coefficient values for perfect source 
eperture correction 
Gaussian ■ — rectangular 

that the taking of the inverse Fourier transform is done by 
expanding the correction aperture spectrum as an even 
power series within the range ±NI2. To a good approxi- 
mation, however, it is only necessary to tal<e the first two 
terms of the series if the amount of correction is low, and 
expressions for the coefficients are obtained as follows: 

Cg = 1 + {7rw)V72 

c, = (-1)''wVl2r^ 

for the rectan^lar aperture of width wjN p.h. 

and "^o = ■" + (^*)^''6 

c, = (-1)'';sVr^ 

for the Gaussian aperture of standard deviation sjN p.h. 

The values w and s can be related to D, the dB of correction 
needed at half sampling frequency, NI2, by the equations 

D=20logjo (1 + 7r^wV24) 
= 20logjo (1 + rr^sV2) 



or 



Note that the simple approximation causes the coefficients 
to fall off as the inverse square, in both cases. 

This result is somewhat surprising as the 'ringing' of a 
siax/x pulse dies away as the inverse. 

2.4. Practical source correction 

2.4.1. Frequency response 

A practical corrector has an aperture of finite width 
which means that frie function of Fig. 3(d), being the 



Fourier transform of the correction aperture, is unbounded. 
It is therefore more helpful to derive the correction coef- 
ficients by way of the transversal-filter response, rather 
than the correction aperture and its transform, because of 
the overlapping of the correction-aperture spectra. Now it 
can be shown that the spectral amplitude response of the 
transversal filter is given by the expression 



/ 



Hm 



+ 2 



r= 1 



c, cosr ( 



where 6 is an angular co-ordinate that varies from -tt to +7t 
across the zero-order spectrum, that is for spatial frequen- 
cies in the range ±NI2 c/p.h. The relationship between d 
and n is thus given by 

e = 2im/N 

If the signal-component amplitudes are restored to their 
correct values the spectral response of the corrector is the 
reciprocal of the scanning-aperture spectrum over the zero- 
order spectral region, that is 

Yy{n)H(e)= 1, |0|<7T 

However, as previously noted, unless the corrector response 
has a cusp at the point = tt and therefore an infinite 
number of terms, it cannot correct perfectly. Thus the 
function Y {n)H(d) cannot equal unity over the whole 
region and the error can be expressed as a deviation from 
unity. The difficulty now arises as to where to allocate 
the error. For example a positive error in one part of the 
spectrum can be offset by a negative error in another part 
without causing undue subjective impairment. The simplest 
yardstick by which to judge the efficacy of correction as a 
function of /, the order of correction, is to allocate the 
error to the higher frequencies by making the overall 
response maximally-flat at zero frequency; this will be 
referred to as the maximally-flat condition. It can be 
achieved by expanding both [{Yyin)]"^ and H{d) as even 
power series and equating coefficients of ascending powers. 
As both power series are infinite this operation produces 
an infinite number of equations containing / + 1 unknowns 
and for the maximally-flat condition the first / + 1 equa- 
tions are satisfied. An example of this process for the 
Gaussian source aperture is given in Appendix 2. When 
the coefficient values are found they can be aibstituted in 
the function H{0) and the overall response Yy{n)H{e) can 
be determined. 

Fig. 7 shows the overall response/frequency charac- 
teristics obtained for values of / up to 4, that is, 4th-order 
correction. A Gaussian source aperture is assumed with 
various values of loss at the half-sampling frequency. As 
can be seen the characteristics never overshoot and the 
effect of adding higher-order correction is to delay the 
onset of the roll-off. There would appear to be a law of 
diminishing returns in that the improvement gained by 
applying first-order correction is more than twice that 
gained by adding second-order correction. Fig. 8 shows the 
corrector response alone for values of / up to 4 and for 10 
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Fig. 7 - Overall response of source aperture and corrector. Gaussian source with corrector adjusted for maximally-flat 

response at zero frequency 
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Fig. 8 - Corrector response adjusted for maximally-flat 

correction of Gaussian source at zero frequency. Source 

loss of 10 dB at half-sampling frequency 

dB of source loss at the half-sampling frequency. Note the 
increasing sharpness of response at the half-sampling fre- 
quency as/ increases. 

2.4.2. Transient response 

Although the respond/frequency characteristic is a 
measure of definition it gives no indication of the transient 
behaviour of the system and in particular tlie degree of 
overshoot which reajlts from correction. The transient 
behaviour can be deduced starting with an appropriate 
signal spectmm and working through the system operations 
in tiie spectral domain, finally transforming back to the 
space/time domain. This process involves the display 
aperture mentioned earlier and the form of the display 
aperture considerably influences the final result, as will be 
seen. 

Fig. 9 shows a computer simulation of the system 
response to a 27 edge; such an edge is defined as the result 
of passing an infinitely sharp edge through a filter having a 
27'-pulse response where T is the period of the highest fre- 
quency that can be displayed by the system. In vertical 
spatial terms MT is equivalent toNI2 and IW corresponds 
to a vertical picture element (PEL). Thus the 23" edge 
under consideration has a width of 2ror 4 PELs. Such an 
edge has been chosen because it contains very little energy 
above the spatial frequency NI2. A source loss of 10 dB at 



the half-sampling frequency and a display aperture having a 
flat response/frequency characteristic cutting sharply at the 
half-sampling frequency have been assumed. The use of 
this display aperture, although physically unrealisable as it 
is of (sin x)lx form which requires negative luminance, is 
nevertheless helpful in indicating the general behaviour of 
the system in response to different degrees of correction 
without introducing unnecessary detail. The correction 
coefficients have been derived according to the criterion 
explained above giving a maximally-flat response at zero 
frequency. 

Because the signal spectral energy is largely confined 
below half the sampling frequency there is virtually no 
difference between the original signal and a perfectly sam- 
pled and displayed version of it as can be seen from Fig. 9. 
If, however, a source loss of 10 dB is introduced and re- 
mains uncorrected the response is appreciably softer. 
Using the '10-90%' criterion, the rise-time is increased by 
some 20%. As can be seen the introduction of first-order 
correction practically restores the rise-time to its former 
value with a 0-6% overshoot. After 4th-order correction 
the rise-time is indistinguishable from the original but an 
overshoot of 0-3% remains. The remaining curve indicates 
the response when the first-order correction is increased to 
the point where the correction at the half-sampling fre- 
quency is equal to that of 4th-order correction (based on 
the criterion of maximally-flat response). This condition 
is indicated in Fig. 8. Thus the over-corrected first-order 
and 4th-order responses could perhaps be said to have the 
same definition for high frequencies. But clearly the over- 
corrected first-order response gives much more overshoot 
than the 4th-order. This exercise illustrates the importance 
of allocating the response error to the higher frequencies. 

Figs. 10 and 11 show the system responses to half- 
sinusoidal edges of widths 3 and 2 PELs respectively. The 
spectmm of the first edge has a zero at the half-sampling 
frequency. The same conditions as in Fig. 9 have been 
assumed. Here, even with perfect sampling and display, 
the response rings because the signal contains significant 
spectral energy beyond the half-sampling frequency. The 
effect is partly caused by aliasing in the sampling operation 
but predominantly by the sharp cut in the assumed display 
response. In Fig. 10, as before, first-order correction makes 
a significant improvement which is little bettered by 4th- 
order correction. For the edge in Fig. 11, which clearly 
defeats the system, first-order correction would probably 
be preferred to 4th-order correction because it reduces the 
overshoots although improving rise-time considerably. 

Lest the foregoing examples seem suspect because of 
the physically unrealisable display aperture. Fig. 12 shows 
the same conditions as Fig. 11 but with a Gaussian display 
aperture having a loss of 10 dB at the half-sampling fre- 
quency. The same general conclusions apply although the 
issue is confused by the presence of the perturbation at the 
sampling frequency. This perturbation is present because 
the display aperture attenuates the sampling frequency by 
only 40 dB instead of completely rejecting it. 

In summary a maximally-flat response at zero fre- 
quency would seem to result in least overshoot on edges 
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Fig. 9 - System response to 2Tedge with various orders of correction. Gaussian source with loss of 10 dB, sharp-cut display 
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Fig. 1 1 - System response to half-sinusoidal edge of width 2 PELs. Other conditions as in Fig. 9 



that are capable of being displayed by the system within 
tile confines of the responses obtainable by finite-order 
correction. Little is gained by increasing correction 

beyond first-order. This degree of correction is also prob- 
ably best for edges which are too sharp for the system and 
which consequently ring when perfectly corrected. 

2.5. The effect of interlace 

The previous discussion has been concerned with a 
system involving sequential sampling. The argument need 
only be slightly modified, however, to deal with interlaced 
sampling. For simple interlaced sampling the spatial 
sampling frequency in each field is halved but the odd 
^ectral orders of the interlaced sampled field are inverted. 
Consequently when the spatial spectra of the two displayed 
fields are combined the odd orders partially cancel to an 
extent which depends on the temporal response of the 
display and the observer. Thus the odd orders are attenu- 
ated and alternate during a complete picture and, taken 
over the integration time of a single picture, they can be 
ignored. The effective spatial sampling frequency over a 
complete picture is therefore equal to that of sequential 
sampling. If, however, aperture correction is performed 
using a transversal filter having access to the samples of 
only one field the spatial sampling frequency appropriate 



to the corrector response corresponds to that of the field. 
Thus the periodicity of the response in the spatial frequency 
domain is halved and is given by 

/ 
//(0) = Cq -t-2 y CyCos2re , 9 = 2Tm/N 

r= 1 

which has a period of NI2 where, as before, a picture (2 
f iel ds) contai ns A^ I ines. As 



/ 



H(0) = c 






c^= 1 



/■= 1 



for unity gain at zero frequency it is clear that the corrector 
offers no correction at the frequency NI2 c/p.h., the effec- 
tive half-sampling frequency for the picture. In fact the 
corrector exerts its maximum effect at the frequency 7V/4 
c/p.h., the mid-band of the effective display response. If 
therefore the device attempts to correct a Gaussian source 
aperture up to this frequency on a maximally-flat basis the 
simultaneous equations which govern the correction coef- 
ficients are obtained, as before, by equating ascending 
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Fig. 12 - As in Fig. 1 1 except Gaussian display aperture with 40 dB rejection of sampling frequency 



powers but halving the value of s, the normalised standard 
deviation. 

If, as before, the overall response is maximally-flat at 
zero frequency the first / + 1 equations are satisfied. Fig. 
13 shows the overall response/frequency characteristics for 
correction up to 4th-order, and for 10 dB source loss at the 
half-sampling frequency. Note that there is now no correc- 
tion at this frequency. Fig. 14 shows the corrector 
re^onse alone. 

Figs. 15, 16 and 17 show the time domain responses 
to 2rand half-sinusoidal edges corresponding to Figs. 9, 10 
and 1 1. The same conditions are assumed but the corrector 
response is now that of Fig. 14. It will be observed in Fig. 
15 that the first-order corrected signal overshoots a little 
more (0-8%) than in Fig. 9. The 4th-order corrected signal 
overshoots even more however and is accompanied by 
further ringing extending away from the edge. This con- 
dition would be likely to be more objectionable. The 
first-order overcorrection is not as extreme as in Fig. 9 as 
might be expected because of the inability of the corrector 
to affect the higher frequencies. For the sharper edges of 
Figs. 16 and 17 the difference between perfect and 4th- 
order correction becomes more marked than in Figs. 10 
and 1 1 whereas the difference between first and 4^h-order 
correction decreases. The overshoots are also worse than 
for the sequential case. 

In summary the law of diminishing returns seems to 



apply more strongly for correction involving interlaced 
sampling than for that involving sequential sampling. 

2.6. The effect of source noise 

As the complexity of the correction increases, more 
and more signals contribute to the output of the corrector. 
These signals inevitably contain noise and assuming the 
noise components are uncorrelated their powers rather 
than amplitudes add together. Thus the signal-to-noise 
ratio at the output of the corrector is less than that at the 
input and the degradation can be calculated as follows. 

If the r.m.s. noise voltage at the input is gj the output 
r.m.s. noise voltage, e^, is given by 

eo'=Co'ei'+2c,2ei^ + 



(NjC""""^ \ 



r= 1 



This noise will be most visible on areas of uniform lumi- 
nance corresponding to a constant signal input. Under 
these conditions the corrector gain is unity and the ratio of 
the input to output signal-to-noise ratios is simply the ratio 
of the noise powers. Thus the corrector noise factor in dB 
is given by 



fPH-118) 



- 10- 




0-1 0-2 0-3 

normalised frequency 



0-4 



0-5 



Fig, 13 - Interlace corrector. Overall response of source 
aperture and corrector. Gaussian source witti loss of lOdB 
at half picture sampling frequency. Corrector adjusted for 
maximally-flat response at zero frequency. Various 

orders of correction 



\ r= 1 / 



F= lOlog 



This relationship is plotted In Fig. 18 for various 
orders of correction and varying degrees of correction. As 
a corrector can be given any noise factor by appropriate 
choice of coefficients the criterion of comparison in Fig. 1 8 
is the maximally-flat condition assuming a Gaussian source 
whereby the coefficient values are determined by the 
degree of correction needed at the half-sampling frequency. 
Although the coefficientvalues fall off rapidly as r increases, 
a higher-order correction allows the values of all coefficients 
to rise for the same degree of correction. Thus for 10 dB 
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Fig. 14 - As Fig. 13 - Corrector response alone 

of correction the value of the main term c^ rises from 1-23 
for a first-order correction to 1-56 for perfect correction. 
As a result the addition of higher order terms does 
noticeably affect the noise factor of the corrector. How- 
ever, the increase in noise factor between third and fourth- 
order correction is much less than that between first and 
second; but fourth-order correction is still much less noisy 
than perfect correction. If a figure of merit is defined as 
the difference in dB between the maximum correction 
offered (as deduced from Fig. 8) and the noise factor for 
the same degree of correction, this quantity gives a more 
meaningful comparison between orders. In fact the quan- 
tity increases with increasing order of correction but once 
again a law of diminishing returns applies. 

Inajmmary.asthecorrection becomes more complex, 
the increase in noise follows roughly the same pattern as 
the increase in correction, so that, other things being equal, 
there is little to choose between various complexities. If, 
however, the noise at the input were always imperceptible 
there could be a clear benefit in introducing more complex 
correction. 



3. Practical tests 

3.1. Introduction 

The foregoing theoretical analysis implies that if the 
source aperture is corrected so that the spectral response is 
maximally-flat at zero frequency there is little to be gained 
by introducing higher-order terms to a first-order correction 
process. It may be, however, that the maximally-flat con- 
dition may not be, subjectively, the best form of correction. 
Accordingly two methods of approach suggest themselves. 
(A) The necessary amount of correction for different orders 
can be predicted according to the maximally-flat condition 
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Fig. 15- Interlaced systBm response to 2T edge. Conditions as in Fig. 9 
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Fig. 16 - Interlaced system re^jonse to half-sinusoidal edge of width 3 PELs. Conditions as in Fig. 9 



knowing the source characteristics, and different orders can 
be compared subjectively. (B) The subjectively optimum 
conditions for each order of correction can be compared. 
Both these methods were tried in an attempt to investigate 
the benefit of adding higher-order correction. 



3.2. Colour considerations 

A preliminary investigation was carried out to test 
whether the results of correcting monochrome pictures 
could be applied to the correction of colour pictures. A 
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Fig. 17 - Interlaced system response to half-sinusoidal edge of width 2 PELs. Conditions as in Fig. 9 

4-tube plumbicon camera was used as the source and a 
digital vertical-aperture corrector was inserted in the feed of 
linear luminance signal to the processing circuits. The 
luminance camera tube was misregistered by an amount 
corresponding to the delay introduced by the corrector so 
that the corrected luminance signal was registered with the 
signals from the other three tubes. Two picture subjects 
were used. The first was a bunch of highly coloured 
flowers containing much fine detail. The second was an 
isolated flower with a thin, near-horizontal leaf against a 
light background. Five different states of correction were 
used, shown in Table 1, and six observers were asked to 
10 -12 -14 16 -18 20 score the pictures, for definition and overshoot separately, 

source loss.dB according to the EBU impairment scale shown in Table 2. 

The colour was then removed from the picture (by shorting 

Fig. 18- Corrector noise factor as a function of source loss together the display guns to avoid cross-colour) and the 

for various orders of correction observers were asked to score the same states of correction. 




TABLE 1 



Correction Characteristics Used in Colour Tests 



Correction 
Condition 




1 
2 
3 
4 



Coefficients 



1-00 

1-23 —155 

1-34 —17 

1-90 —60 

2-70 —85 




•f-04 


-1-15 





dB loss at half-sampling frequency 
Field Picture 





8 

16 

24 

>40 





32 

64 

96 

>160 
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Fig. 19- Subjective grading of correction conditions for colour and monochrome pictures 



TABLE 2 
EBU Impairment Scale 

Grade Degree of Impairment 

1 Imperceptible 

2 Just perceptible 

3 Definitely perceptible but not disturbing 

4 Somewhat objectionable 

5 Definitely objectionable 

6 Unusable 

The results, averaged over observers, are shown in 
Fig. 19. The dashed lines indicate the standard deviation 
range. Table 1 shows that correction condition 4 corres- 
ponds to greater than 40 dB source loss at the half-field 
sampling frequency. As this is much greater than the loss 
of even the worst plumbicon tubes this condition would 



be expected to overcorrect the pictures grossly. In fact the 
observers rated this condition as giving the best improve- 
ment in definition whilst giving the worst overshoot. The 
other conditions as might be expected, were rated in 
between this limit and that of giving no improvement. As 
can be seen there is a high inverse correlation between defi- 
nition and overshoot and a high correlation between reailts 
for colour and monochrome pictures. When calculated for 
individual observers the self-correlation coefficient for 
colour and monochrome grading never fell below 0-9. It 
was thus concluded that this was sufficient justification for 
applying the results of monochrome picture tests to colour 
pictures. 

3.3. The measurement of source aperture 

The correction prediction method mentioned above 
requires a knowledge of the source aperture profile. It was 
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decided to use a high-quality slide scanner as the source so 
that the source aperture, being ihe CRT spot profile, was, 
to some extent, under the control of the experimenter and 
would be relatively easy to measure. In fact there are many 
pitfalls in carying out this meaajrement and it. proved" 
niore difficult than anticipated. 

There are three common ways of measuring a C.R.T. 
^ot profile apart from direct measurement by microscope. 

(a) The horizontal frequency method. In this method 
the spot scans a grating of known pitch and the light 
is collected by a photodetector whose output gives a 
measure of the response/frequency characteristic of 
the aperture at the grating frequency. By using dif- 
ferent grating pitches the spectrum of the aperture 
can be built up and the Fourier transform gives the 
spot profile. If the grating is not sinusoidal care 
must be exercised in interpreting the photodetector 
output The disadvantage of this method is that the 
spot must scan sufficiently slowly to eliminate errors 
caused by afterglow response in the C.R.T. phosphor. 
If, at normal scanning speeds, the afterglow correction 
is at all suspect, this method cannot be used. 

(b) The knife-edge method. In this method the spot is 
scanned slowly across a knife-edge and the output 
of a photodetector collecting the transmitted light is 
proportional to the integral of the spot profile. Thus 
by differentiation the profile is obtained. The diffi- 
culty with this method is the maintenance of mechani- 
cal stability needed to produce repeatable results. 

(c) The vertical frequency method. Here a spot scanning 
horizontally at normal television line frequencies is 
slowly scanned vertically across a horizontal grating 
and Ihe photodetector output gives a spot measure- 
ment of the response/frequency characteristic as in 



(a). The difference between this method and that of 
(a) is simply the speed of scanning and the benefit of 
integration conferred by the horizontal scanning. 

All these methods were tried and method (c) was 
found to be the most reliable. The vertical scanning was 
performed manually by operating the shift control on the 
scanner. The normal horizontal and vertical scans were left 
in operation and by strobing out a particular line of the 
photodetector output waveform the slow change in magni- 
tude of the line signal could be observed as the shift was 
operated. The peak-to-peak magnitude of the signal gave a 
measure of the response at the grating frequency. Several 
grating frequencies were used (the gratings on Test Card C 
were found convenient) whose pitches were accurately 
measured. The results could be expressed in cycles/mm on 
the tube face knowing the magnification ratio of the optics 
associated with the scanner. 

Fig. 20 shows the results obtained for various posi- 
tions of the electronic focus control. The area of uncer- 
tainty accompanying response values near unity comes 
about because a square-wave grating was used. With such a 
grating a correction factor of 4/7r must be applied to peak- 
to-peak measurements if the waveform is square but not if 
it is sinusoidal. The curves of Fig. 20 were subjected to 
best-fit analyses assuming Gaussian spectra and values for 
the standard deviation of the spot profile in /im were de- 
duced shown in Table 3. 

Referring to Section 2.3 the amount of correction 
required depends on the quantity s, the ratio of the spot 
size to the scanning pitch. The scanning pitch was found 
by measuring the number of field lines corresponding to a 
certain measured distance on a slide and converting this to 
a measurement on the tube face. The values of s appro- 
priate to the field scanning pitch thus deduced are shown 
in Table 3. This quantity can be used in the simultaneous 
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TABLES 
Source Characteristics 



Source 
Condition 


Spot standard 
deviation (/xm) 


Sf 


dB loss at half-sampling 
frequency 




Predicted correction 
coefficients 










field 


picture 


^0 


^1 ^2 


S 


4 


112 


■66 


18-6 


74-4 


1-436 
1-693 
1-859 


-218 

—390 -043 

— 513 -092 


—008 


5 


67 


•40 


6-85 


27^4 


1-163 
1-224 
1-255 


-082 

— 122 -010 

— 145 -019 


-002 


6 


55 


•325 


4-5 


18 


1-104 
1-138 
1-156 


-052 

—075 -006 

—088 •Oil 


-001 


7 


94 


•555 


13^2 


52^8 


V308 
V458 
1-544 


-•154 

-•254 025 

—319 •OSI 


—004 


8 


156 


•925 


36^5 


146 


V849 
2-603 
3-254 


—425 

—927 -126 
-1-415 -321 


—033 



equations to predict the correction coefficients because the 
transversal-filter delay element corresponds to the pitch of 
a field line. Also shown in Table 3 are the source losses at 
the half-field sampling and half-picture sampling frequen- 
cies. As the worst plumbicon tubes have.a loss of 20 dB at 
the half-picture sampling frequency' it can be seen that 
only positions 5 and 6 correspond to probable source 
characteristics. 

3.4. Subjective assessments of predicted correction 

The coefficients necessary to correct the various 
source apertures according to the criterion of maximally- 
flat response are shown in Table 3. As the coefficients in 
the digital aperture corrector were quantised to 7-bit 
accuracy it would have been possible to apply 3rd-order 
correction for only conditions 4 and 8. Consequently it 
was decided to concentrate on differences between 1st and 
2nd-order correction. Source conditions 4 and 5 were used 
and for each condition observers were asked to judge the 
difference between zero (no correction) and first-order, and 
between first and 2nd-order correction, according to the 
EBU Comparative Scale shown in Table 4. Three slides 
were used as source material; the first was a Test Card, the 
second a scene having random fine detail and a black out- 
line (boat), and the third contained more ordered detail 
(Maida Vale Studio). 

The results are shown in histogram form in Fig. 21. 
The hatched area represents the standard deviation and a 
single line indicates a unanimous score. For source con- 
dition 5 only Test Card was judged to be significantly 
improved by adding 2nd-order correction and only by a 



TABLE 4 



EBU Comparative Scale 



Grade 


Degree of Comparison 


+3 


Much better than 


+2 


Better than 


+ 1 


Slightly better than 





Same as 


-1 


Slightly worse than 


-2 


Worse than 


-3 


Much worse than 



grade. For source condition 4, Test Card and Maida Vale 
were judged to be improved by a grade. The general 
comment was that insufficient correction had consistently 
been applied. 

3.5. Subjective assessment of optimally corrected pic- 
tures 

Following the results of the above test it was con- 
cluded that the criterion of maximally-flat response was 
not subjectively the most pleasing. A further test was 
therefore devised in which observers were first asked to 
indicate the degree of first-order correction they preferred. 
Two source conditions and two slides were used. This 
degree of correction was then held as a reference and the 
observers were asked to judge between this reference and 
various degrees of second-order correction using the EBU 
comparative scale. The amount of second-order correction 
was varied within a range which included the degree of 
first-order correction for low frequencies. 
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Fig. 21 - Results of correction prediction subjective tests 

— 1 Compa-ison between no correction and first-order correction 

1 — 2 Comparison between first order and second-order correction 
Note: Tlie sign of the comparative grade is positive if the last 

named condition is judged better than the first condition 



The results of the first part of the test are shown in 
Table 5. The figures for losses at the half-field sampling 
frequency are those for the source, calculated from the 
coefficients according to the maximally-flat criterion. 
Comparing these figures with those of Table 3 calculated 
from actual source meaajrements it will be seen that there 
was always a tendency to over-correct although less so with 
the Test Card. Fig. 22 shows the response/frequency 
characteristic resulting from the coefficient values of Table 
5 in conjunction with the source spectra calculated from 
source measurements. 

This tendency to over-correct is not entirely under- 
stood but could have been caused in part by the display 
aperture which unfortunately was not measured. A display 
aperture spectrum with excessive attenuation at high 
spatial frequencies would require compensation in the 
signal fed to the display. 

The results of the second part of the test are shown 
in Fig. 23 in which comparative grade (with respect to the 
standard) is plotted as a function of the two independent 
second-order coefficients, Cj and c^. It will be seen that 
for the great majority of conditions, second-order correc- 
tion was judged inferior to first-order correction. 
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Fig. 22 - Preferred overall response for two scenes and two 
source conditions 



TABLE 5 



Sou rce 


Scene 


Correction 


dB loss at half-sampling frequency 


Condition 




Coefficients 


field picture 


6 


Boat 


1-43 -215 


18-4 


73-6 


5 


Boat 


1-46 —23 


19-7 


78-8 


6 


Test Card 


1-19 —095 


8-1 


32-4 


5 


Test Card 


1-26 —13 


11-1 


44-4 
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Fig. 23 - Subjective assessment of second-order correction by comparison with first-order preferred correction, (a), (b) and 
(c) represent different source conditions and picture materiai. Comparative grade ispiotted on a three-dimensional graph as 
a function of the two independent second-order variables. The plotted points give the values of c j and c^ appropriate to the 
second-order correction conditions and the numbers beside the points are values of the subjective grading of the conditions. 
Note the grading is positive if the condition is judged better than the standard. The dotted line indicates the relationship 
between c^ and c^ when the response/frequency characteristic of the overall system has the same behaviour near the 



frequency origin as the standard condition 



4. Conclusions 



It would appear from both theoretical and practical 
considerations that there is little to be gained from increas- 
ing the complexity of aperture correction beyond first- 
order. A relatively crude analysis of perfect correction 
reveals that the correction coefficient values decrease as the 
inverse square of the distance from the wanted signal. An 
analysis of practical correction indicates that the rate of 
decrease is even faster. Consequently the higher-order 
terms rapidly become insignificant. This law of diminish- 
ing returns is confirmed by the theoretical time-domain 
response to various edges and especially in the case of 
interlaced-sampling where higher order correction may 
a;tually impair the response by introducing ringing before 
and after edges. Subjective tests indicate that there is a 
slight preference for second-order correction over first- 
order if the correction is of such magnitude that the fre- 
quency response is maximally-flat at zero frequency, and 
does not overshoot However, first-order over-correction is 
preferred to either of these correction conditions. 

These subjective results are based on tests using a 
corrector operating on the signals related to a single field. 
If, in an interlaced system, it were possible to operate on 
the signals of both fields there might be a slight gain in 
introducing second-order correction but this would only 
apply for stationary scenes. Even for these Kienes, using 
ajch refined correction, there would probably be difficulty 
with 'twittering' on horizontal edges. This is purely a 
defect of conventional display devices which rely on the 



fact that in an interlaced system the picture content does 
not exploit the full spatial bandwidth of the two combined 
fields. 

For moving scenes, at least a picture store would be 
required in order to interpolate points in the interlaced 
field and the present cost of thjs process would probably be 
deemed excessive in relation to the benefit gained. 
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Appendix 1 

The Derivation of Correction Coefficients for Perfect Source Correction 

Following the method outlined in Section 2.3 it is assumed that the correction-aperture spectrum can be expressed 
as an even power series in the ^atial frequency variable, that is 

F{v)=N-^['\ + {u/p^)^ + {v/v^f + ], \v\<Nn 

= otherwise 

where A'^ is the spatial scanning frequency. The Fourier transform can be written 






fe= 



where/, J. results from the term in v'^'^. In general 



-N/2 
Now the identity 



f2k^^ = ^~''lt^2k^^ Where 

iV/2 



V t fc d^ /2sin A 

J/exp-iy.d.= j^^|-y-j 



is derivable from a standard general form. Thus 



I^^ty)- J ^^'^exp2.j^d.-A^^--^^— J~^I 

-Nl2 ^ ' 

Jn^ d^^ /sin z\ 
\2/ dz^^ \ z I 



€ 



Af(-ir 



sinz 2A;cosz sinz ,,/t,„, ., JlHi. 
+ _^ _ 2k{2k - 1)-^- ... + (-1)^^(2^)! -^-^ 

z z^ z z 



where z = Af7rF 
The coefficients are points on the correction aperture spaced at internals of MN, that is at the points r/N where r is an integer. 

oo oo 

k= k= 
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and 



Forr= 



Thus 



Ldr/N) = N 



■■N 






stnrrr IkcosTrr slmrr , sin rrr 
+ :; 2k(2k- 1) - . ..+ (-1)*(2yt)! r-: 



(rrr)' 



2A: 2A:(2A:-1)(2A:-2) 



(irr)^ (wf 



+ ....+ 



(-1)*-^2A;)! 



(7Tr) 



2 A: 



(■nr? 



, r¥=Q.k¥=0 



',.(0, = A-,-„*(|) „.^|,_i!.|_....J 



/o(0) = 7V 



/2(0) = 7V(-1) 



V / '^ V 



-2 4.3 z^ 

— + 

3! 5! 



) 



2/3! 



/^(0)=7V(-1)M — J lim 



.2A: 

/,,(0)=7V|-™ 



^4! 6.5.4.3.Z2 \ _ /tvV 4! 



2 / 2A: + 1 



Therefore 



,^-^ 1 f ^ \ 



(-ir 



°° T/ \ 2fe k 

El N \ Y^ (2 
\2p^, I L^ [2k + 1 



2k 2A;(2A;-1)(2A:-2) 



(2^)! (-1)'^- 



- 2q)\ (w?'^ 



For the Gaussian aperture of standard deviation a 

F[.v) = TV- • exp 27T^ a^ j^^ li'K AV2 



But 



so that 



otherwise. 



exp 27r^CT^i^^ = 1 + 27r^a^i^^ + (27r^a^i^^)^/2! + 



■Ik 



v^^ -(2TT^a^f/k\ 



(PH-118) 



- 20 



Thus 



oo / V 2 ^ OQ 

2k + 1 A:! I 2/ Z_^(2fc+1 

it=0 ^' A;=0 



)fc! 



and 



c, = (-ir 



,2k ^-^ ^2k)\ (-1 )'?-■' 



fc= 1 






L <7= 1 



2q)\ (■nr?'^ 



where z = ■naN/\/2 = n x/\/2 

For the rectangular aperture of width w/N 

■nwv/N 



F{u) = TV 



I / irwv/N \ 
I sin nwv/N I 



\i>\<N/2 



sin TTWv/NJ 
= otherwise 

Now it can be shown that the expansion of x/sinx has the form 
x/sinx = 1 + A.x^ + A.x"* + 



H<2 



'2 4 



where 



Aj = 1/3! 



A^ = 1/(3!)^ - 1/5! 



A, = 1/(3!)' - 2/3!5! + 1/7! etc. 

6 



The ratio of successive coefficients tends to tt 



Thus 

and 
Thus 

and 



F{v) = Af~ M 1 + Aj (ttwv/N)^ + A^ (tnvv/N)'* + 



p;^=(/V/7rwpVA,fc 

OO 

/t= 



2/t 






(2/1:)! (-I)*? 



-1 



-2^)! {Trr)^'^ 



where z = 7rw/2. 

The values of z can be related to D, the dB of source loss at the frequency N/2 by the equation 
£»=20log,JF(7V/2)/F(0)] 
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(la) 



(lb) 



(2a) 



(2b) 



For the Gaussian source 

£)=20logjQ exp 27r^a^jVV4 
= 8-68.... Ti^s^n 

= 8-68 z^ (3) 

For the rectangular source 

D = 20 loQjQ [ (7rw/2)/sin(rrw/2)] 
= 20log,Q (z/sinz) (4) 

By eliminatingz in Equations (1) and (3) or (2) and (4) the coefficient values can be deduced as functions of source loss. 

Appendix 2 
The Derivation of Coefficient Values for Finite-order Correction 

For finite-order correction on a maximally-flat basis ascending powers of [Ky («)]"' and //(0) are equated. 
Assuming a Gaussian source we have 



y 
For /th order correction 



[rj«)]-' =exp27r2a2n^ 



/ 



Hm = Cg + 2y c^cosrd 
r= 1 

where d = lirn/N 

Thus equating power series 



■ Cg -I- 2Cj (1 - 0V2I -t- e-'M! -H . . . ) -H 2C2(1 - 2^0^21 -f 2''0''/4l + . . . ) 



H- 



which leads to the simultaneous equations 

Cg/2-)-C, +^2 -H -I-C, = y2 

c, -^2^0^ +....-H/^c, = -% (21/11)2- MaAf)^ 

c, -H2^^C2 + ... + /^''c/ = y2(-1)''[(2^)!/?!] 2-f(aN)'^f 

etc. 

Thus we have an infinite number of equations with only/ + 1 unknowns. On a maximally-flat basis the first/ + 1 equations 
are satisfied giving values of the coefficients c^ . . . Cj as a function of the parameter aN (= s, the normalised standard 
deviation). 
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